Wavepacket-current < — ► supercurrent conversion in closed N/S/N circuits 
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The rise and decay in space and time of a supercurrent 
in a superconductor between two normal current leads con- 
nected to a reservoir is computed with the help of quasiparti- 
cle wavepackets which suffer Andreev scattering in the inter- 
faces between the normal current leads and the superconduc- 
tor. Stills from a computer movie visualize electron *-+ hole 
scattering and the creation and destruction of the supercur- 
rent. 

PACS: 74.25.Fy, 74.80.Dm, 74.80.Fp 



Andreev scattering (AS) of electrons into holes and 
vice versa by spatial variations of the superconducting 
pair potential H, in competition and cooperation with 
conventional scattering, determines the electronic struc- 
ture and transport properties of inhomogeneous super- 
conductors. The Tomasch effect in tunnel junctions 0, 
Josephson currents p,[4|j5|Jfj, 7|fl . ex cess currents, and sub- 
harmonic gap structures [p |l(||lT| , |l2| ] in superconducting 
(S')-normal conducting (A)-supcrconducting junctions, 
as well as the transfer of half of the Magnus force to the 
core electrons of a moving vortex line |13|| are due to AS. 
There is AS in He 3 , too pi . A wealth of AS phenomena 
has been discussed recently in [tta and |16[ . 

While the conversion of a normal current into a su- 
percurrent by electron — > hole scattering in the interface 
between an N and an S region of semi-infinite lengths 
has been described before ]l7| , |l8| ,|9[ , it is the purpose 
of this paper to analyze the normal-current < — ► super- 
current conversion processes in a superconducting layer 
of finite length L z between two normal current leads. 
These normal leads are connected to a reservoir ("bat- 
tery") which acts as the current source in the closed cir- 
cuit. The extensions of the N and S regions in x- and 



y-directions are L x and L y . The metal- vacuum bound- 
aries are treated as rigid walls. By varying L x and L y 
one can vary the dimensionality of the system. The 
much investigated quasi-two-dimensional (Q2D) super- 
conducting/semiconducting heterojunctions |19Lp0f and 
superconducting quantum dots in Q2D channels are sys- 
tems for which our analysis may prove useful. 

AS and the associated formation and destruction of 
Cooper pairs and supercurrents can be calculated from 
the Time-dependent Bogoliubov-de Gennes Equations 
(TdBdGE) 0,|l|,|l],§. They describe the evolution 
of the spinor quasiparticle (q.p.) wavefunction with 
the electron component u n (r, t) and the hole component 
v n (r, t) under the influence of scalar and vector potentials 
V(r, t) and A(r, t) in the single-electron Hamiltonian 
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Here, the matrix Hamiltonian TC(r,t) has Ho(r,t) and 
— iJg(r,t) in the diagonal, and the pair potential A(r, t) 
and its complex-conjugate in the off-diagonal. The chem- 
ical potential /i in Ho is that of the reservoir. We neglect 
all influences of entropy production associated with cur- 
rent flow on the chemical potential, because the number 
of degrees of freedom of the reservoir is assumed to be 
much larger than that of the normal leads and the super- 
conductor. Thus, /i is constant in space and time [gjfl . 

AS is a many-body process. For its analysis it is con- 
venient to consider a non-equilibrium configuration |TJ CT ) 
of the many-body system where one quasiparticle state 
(la), characterized by a tripel I of quantum numbers and 
spin <7, is definitely occupied and all other q.p. states 
(na) are occupied according to the equilibrium Fermi dis- 
tribution function f n . All interactions that might affect 
the spin are neglected. Then, it has been shown jl7 21 
with the help of the TdBdGE (0) that the expectation 
values (T la \p{r,t)\T la ) and (T la \j{r,t)\T la ) of the many- 
body charge- and current-density operators satisfy the 
relation 



d_ 

di 



{T la \p(r,t)\T l(T } + div(T l(T \i(r,t)\T h 



= -4~Im[A*(r,t) Ul (r,t)vt(r,t)} (1 - /,) + divj s j. 
n 

The electron and hole wavefunctions ui and vi satisfy 
eq. (m), and j s ; is the supercurrent density induced by 
the momentum transfer from the q.p. in (la) to the su- 
perconducting condensate by AS. The requirement that 
charge is conserved in the many-body system results in 
the source equation for the supercurrent density 

divj si = 4^Im[A*(r,t)«,(r,t)t;f(r,t)](l - /,). (2) 

A shifted Fermi sphere (or its equivalent in Q2D and 
Q1D conductors) represents the current-carrying many- 
body configurations in the two parts of the normal cur- 
rent leads that are parallel to the z-axis and connected to 
the superconductor. (We do not consider the parts of the 
normal current leads bent towards the reservoir.) In this 



non-equilibrium distribution of electrons above the Fermi 
surface in states with positive momentum in z-direction 
and unoccupied states with negative ^-momentum below 
the Fermi surface the force of the electric field from the 
battery is balanced by the frictional forces from electron- 
phonon interaction. In the following we try to obtain the 
details of normal current < — ► supercurrent conversion 
by studying the motion of electron and hole wavepackets 
that are part of the shifted Fermi sphere. The quasipar- 
ticles in this normal-state configuration are uncorrelated. 
Thus, the total current in the closed circuit is the sum of 
the currents from the individual wavepackets. 

The z-momenta of the electrons (+) and holes (— ) 
in the shifted Fermi sphere are Tik , with fc ± (i?) = 
[k 2 zF ± E2m/h 2 } 1 / 2 . Here, k zF = [k% - {n x n/L x ) 2 - 
(jiy'K j ' Ly) 2 ] 1 / 2 is the z-component of the Fermi wavenum- 
ber &p = pm/x/ft 2 ] 1 / 2 , and (n x ir/L x ) and (n y ir/L y ), n x . y 
integers, are the wavenumbers of the standing waves be- 
tween the rigid walls that limit the metals in x- and y- 
directions. The energy E of both electrons and holes is 
positive and measured relative to the surface of the un- 
shifted Fermi sphere at the chemical potential fj,. For nor- 
mal current densities below the critical current densities 
of conventional superconductors all E are less than the 
modulus A of the pair potential A(z) w A-8(z)0(L z — z), 
where @(z) is the Heavyside function. 

The construction of the representative wavepacket en- 
semble, whose motion and conversion into a supercur- 
rent models the elementary process of current flow in 
the N/S/N circuit, starts with a normalized electron 
wavepacket, localized around zq < in the normal cur- 
rent lead to the left of the superconductor at time t = 0. 
In the center of the wavepacket the energy is E[, We 
choose a Gaussian spectral function 



D(E) 
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the position- uncertainty parameter a z <C \zq\ is chosen so 
large that the related energy spread of the wavepacket, 
SE = H k z p/ma z , is less than A — E\. Solutions of eq. 
(|j), where V and A are neglected and A(r, t) is approxi- 
mated by the real A(z) (thereby neglecting repercussions 
of the q.p. -induced supercurrent on the q.p. and on itself, 
assuming sufficiently small L x and L y ) are multiplied by 
the spectral function D(E), integrated over all energies, 
and matched at the left N/S interface, i.e. z = 0, in the 
usual Andreev approximation, neglecting terms of the or- 
der of A//x outside the exponentials. In so doing, energy- 
dependent functions are Taylor expanded around Ei up 
to first order in (E — Ei). This affects especially k (E) 
and the amplitude of the Andreev-reflection probability 
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Ef] x l 2 is the time for one electron 



where t\ = h[A 2 

— ► hole-scattering event and the associated formation of 

a Cooper pair, see eqs. (||,|]). The resulting electron and 



hole wavepackets Ujvx(r, £) and ujvL(r,£) in the left nor- 
mal current lead, z < 0, and the exponentially decaying 
solutions Msi(r, t) and vsl(t, t) in z > 0, that contribute 
to the source equation (0) essentially in the left half of 
the superconductor, turn out to be 



vnl = wi-y{E l )e ik ^ z e-^ 0+ ^/ v r>+ v t( t -n )? 1^1 

a ik zF z —Kiz -[zo+v+t] 2 /2a1 
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Ef) 1 ' 2 /Tiv z f = 1/tiv z f, with v z p = fik z pjm. [For the 
sake of simplicity the complex wavenumbers in usl and 
vsl have not been Taylor expanded in (E—Ei) but rather 
taken at Ei right away] 

Identifying the wavefunctions usl and v$l from eqs. 
(5) and (o) with the m and vi of the source equation 
(2) and integrating that equation from z = to z yields 
the density of the supercurrent in z-direction induced in 
the left half of the superconductor by AS of the electron 
wavepacket unl into the hole wavepacket vml as 



3 S LL=e z (2ev zF )\w l \ 2 [l~e- 2K ' z ] 

xe -{[zo+v+(t-T l /2)] 2 + (T l v+/2) 2 }/a 2 ( l _ ^y /^ 



Here we have assumed that L z 3> 1/kj. In the oppo- 
site case one would have to add a second source term on 
the r.h.s. of eq. (0) which would contain the contribu- 
tion from the solution usR(T,t),vsR(r,t) of eq. (|l|) for 
< z < L z that matches to the current-carrying q.p. 
wavepackets in the right normal current lead at z = L z 
and decays exponentially with increasing distance from 
L z . This solution determines the supercurrent in the 
right half of the superconductor essentially, and its en- 
ergy and amplitude are obtained as E% and wi from the 
requirement that the two supercurrent densities, com- 
puted from uslv* sl and usrv SRi join smoothly at all 
times somewhere within the superconductor. (Since the 
many-body configurations in the left and right normal 
current leads are uncorrelated, only the current densities, 
not the wavefunctions, must join smoothly.) Because of 
the symmetry of the problem the matching point turns 
out to be L z /2, and we obtain 

U S R = Wie ik ^z e -^(L z -z) e -[z +v+(t-n)] 2 /2a 2 ^ (g) 

vsr = w n (E l )- 1 e^ FZ e-^ L '- z ^e-^ Z0+v ^ t ^/ 2a2 . (9) 

The supercurrent density j s l,R that results from integrat- 
ing eq. (0) from L z to z, with usrv sr in the place of 
uiV{ , has the same form as j s i,L of eq. (M) except that 



exp[— 2kiz] is being replaced by exp[— 2m(L z — z)]. Fi- 
nally, the wavepacket solutions unr(j, t), vnr(y, t) of eq. 
(|l|) in the right normal current lead, z > L z , that match 
to the «sfl( r i t), vsr(y, t) at the right interface in z = L z , 
become 



-[z +L z -z+v+ (t-r,)] 2 /2al 



VNR = Wll(El) V 



! e * 



xe -[za + (v+ /v-)(z-Lz)+v+t] 2 /2a 2 ^ 



(10) 

(11) 



Comparison of the ujvl, vnl with the ujv_r, «at,r shows 
that the center of the electron wavepacket wjvl, propa- 
gating to the right with velocity w ; + in the left normal 
current lead, and the center of the hole wavepacket v^r, 
propagating to the left with velocity vf in the right nor- 
mal current lead, hit the left and right interfaces at z = 
and z = L z at the same time to = —zq/v^, while the hole 
wavepacket vnl, propagating to the left in z < with 
vf, and the electron wavepacket unr, propagating to the 
right in z > L z , are retarded by the time 77 with respect 
to the incident wavepackets. This means: electron — ► 
hole scattering in the left interface combines the electron 
with spin up above the Fermi surface with another spin- 
down electron below the Fermi surface to form a Cooper 
pair at the chemical potential, leaving the hole with spin 
up in the left normal region, transferring the momentum 
2hk Z F to the condensate, and simultaneously in the right 
interface a Cooper pair decays by hole — > electron scat- 
tering into the incident hole and the reflected electron 
of total momentum 2fik Z F- The supercurrent density j s ; 
spreads instantaneously in the superconductor. Its spa- 
tial maximum is at z = L z /2, and its maximum in time 
occurs at t = to + n/2. From the equations for j s ; , usl,r 
and vsl,r one sees that the quasiparticle current density 



and when superconductivity has broken down the uncor- 
rected normal-state configuration reigns everywhere in 
the circuit. If the single S layer is replaced by a meso- 
scopic SNS junction, the many-body configuration in the 
central N layer is a phase-coherent one and thus differ- 
ent from the uncorrelated configurations in the normal 
current leads. In an N/SNS/N circuit the SNS junc- 
tion acts as a gapless superconductor [|4J. It can carry 
a dissipation-free Josephson current through the central 
N layer via phasecoherent q.p. states with \E\ < A and 
\E\ > A Ij. This current converts as a whole into the to- 
tal supercurrent of the S layers, and vice versa, whereas, 
according to eq. (0), each uncorrelated q.p. from the nor- 
mal current leads individually induces its proper fraction 
of the total supercurrent. If the total current density ex- 
ceeds the critical Josephson-current density at a Fermi- 
sphere shift of hq c sNS ~ h/d, where d is the length of 
the central N layer M , a voltage drop appears across the 
SNS junction. There are different models @Jpjl|,||,|| 
for SNS junctions with voltage drops. They differ with 
respect to the implicit assumptions about the rate and 
energy range of q.p. creation in the central N layer by 
supercurrent — ► quasiparticle-current conversion. The 
question of how this rate and range depend upon the 
weakening of phase coherence in the SNS junction by en- 
ergy exchange between quasiparticles and electric fields, 
phonons, and thermal fluctuations like Nyquist-Johnson 
noise is presently investigated. 



Jqpi 



-Re 



h h 

U, —VUi — Vi — VVi 

i i 



(1 - //) 



changes into the supercurrent density j s ; and vice versa 
within a distance 1/ Ki = v z ft~i from the interfaces. Fig- 
ure 1 shows stills from a computer movie of the de- 
scribed Andreev-scattering and wavepacket-current < — ► 
supercurrent-conversion processes. 

By summing the current contributions from the 
steady flow of all quasiparticle wavepackets (7c), I = 
(Ei, n Xl n y ), that corresponds to the shifted Fermi sphere 
in the normal current leads, one obtains the total nor- 
mal current and the total supercurrent within the closed 
N/S/N circuit. 

The supercurrent, carried by the condensate in the S 
layer, involves only states with \E n \ > A. It contin- 
ues the current from the electrons and holes with ener- 
gies Ei < A in the normal current leads. If the total 
current density exceeds its critical value, i.e. if the cen- 
ter of the Fermi sphere in the normal current leads is 
shifted by more than Tiq c s — Am/Tikp, depairing sets in, 
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FIG. 1. Propagation and Andreev scattering of the prob- 
ability densities \u\ 2 and \v\ 2 of spin-up electron and hole 
wavepackets (solid lines) and the induced supercurrrent den- 
sity js (dashed line). The energy in the center of the 
wavepackets is Ei — A/2 = 0.15 meV, the spatial wavepacket 
spread is a z — 5 /j,m and k Z F = 0.9/cf = 2.06 nni -1 . The 
time elapsed between the first and the last picture is 160 
ps. The time elapsed between the center of the incident 
left electron (right hole) wavepacket hitting the interface 
at z = 0, (z — L z — 10 (im) and the center of the re- 
flected left hole (right electron) wavepacket leaving the su- 
perconductor is n = h/[A 2 - £?] 1/2 = 2.53 ps. Note 
that current flow in a closed N/S/N circuit always in- 
volves electron — > hole scattering in one N/S interface and 
hole — > electron scattering in the other interface. The An- 
dreev-refiected wavepackets may be considered as supercur- 
rent-transmitted wavepackets as well. More stills from other 
computer movies on electron — » hole and electron — ► electron 
scattering for various Ei and k Z F in one interface between 
semi-i nfinite iV and S layers can be viewed in the Internet 



under http: / ' /theorie.physik. uni-wuerzburg. de/TPl/kuemmel/ 
profile.html. 
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